Effect of the vortex core on the magnetic field in hard superconductors 
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Using approximate analytical and new numerical solutions of the conventional Ginzburg-Landau 
equations we calculate the small angle neutron scattering cross-section and the variance of the field 
distribution as measured by muon-spin rotation for superconductors with large Ginzburg-Landau 
parameter k. Our results prove that a proper account of the finite size of the vortex core is important, 
even at relatively low fields. This finding provides a natural explanation for the recently observed 
field dependence of the CeRu2 form factor and of the YBa2Cu306.g5 penetration depth. 
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The study of the vortex state in high temperature 
and heavy fermion superconductors is presently a sub- 
ject of intense investigation. Numerous publications are 
devoted to the measurement of the magnetic penetration 
length A since this is one way to probe the nature of the 
low energy excitations and the symmetry of the pairing 
state. Among the possible experimental techniques avail- 
able to investigate the vortex lattice, small angle neutron 
scattering (SANS) and muon spin rotation (/iSR) exper- 
iments are unique since they directly probe the bulk of 
the material and allow to determine not only the field 
and temperature dependence of A but also its value at 
low temperature, see the recent Refs. |l|,|,||||,§. To 
extract quantitative information from SANS and /iSR 
measurements, a detailed theory of the magnetic field 
inside the superconductor is needed, going beyond the 
London model which treats the vortex cores as mathe- 
matical singularities. The finite core size was considered 
in Refs. @,|,|. 

In this paper we compute the Fourier components of 
the magnetic field in a type-II superconductor contain- 
ing an ideal vortex lattice. We disregard pinning pC| ] 
and vortex "phases" such as the glassy or liquid states 
. When accounting for the finite size of the vortex 
cores within the Ginzburg-Landau (GL) theory we find 
an unexpected large reduction of all Fourier components 
down to very low inductions B. Although our results are 
based on the conventional GL theory, they still are of rele- 
vance for the analysis of unconventional superconductors 
such as high Tc superconductors and heavy fermion su- 
perconductors. For example, in recent reports ]r^ , p^ the 
effect of the finite size of the vortex core is described as 
if these compounds were conventional superconductors. 

We define an orthogonal reference frame {x,y,z), with 
the external magnetic field Bext applied along the z axis 



chosen along one of the three main axes a, b and c of the 
penetration-length tensor such that the vortices are also 
along z. For superconductors with large GL parameter 
K = A/f ^ 1 (f is the coherence length) at not too large 
fields Bcxt ^ {Bc2 is the upper critical field) we may 
approximate the vortex fields by the London model. The 
London field B(r) caused by straight vortices located at 
sites satisfies 15|T^ 



B(r) -I- curl[AcurlB(r)] = 



(1) 



Here $o = 2.07 x 10~^^ Tm^ is the quantum of flux, 
the sum is over the vortices, i5(r) is the two-dimensional 
delta function, and z is the unit vector along the vortex 
cores. The eigenvalues of the tensor A are expressed in 
terms of penetration lengths: A^ = A^, A;, = A^, and Ac 
~ Ag . Here Aq , Af, and Ac are the penetration lengths for 
currents flowing along the a, h and c axes, respectively. 

When the vortices form a regular lattice it is con- 
venient to introduce the Fourier components B(G) = 
/ B(r) exp(— iG • r)cPr/S of the periodic magnetic field 
B(i") = J2g ^(G) exp(?G • r), where G are the vectors of 
the reciprocal lattice and S the surface of the vortex lat- 
tice unit cell. The London equation is then easily solved 
for the cases of main interest, namely, Bext parallel to 
either a, b or c. For these three geometries one finds 



BAG) = ^ 



1 



S 1 



^xGy 



AyGl 



(2) 



and Bx{G) = By{G) — 0. Therefore, as expected, there 
is no transverse field component. Equation (^) means, 
for example, that if Bext || c we write this equation with 
X ^ a, y = b, and z = c. In this way we recover the result 
of Ref. for a uniaxial superconductor, in which two 
penetration lengths are equal. 
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Equation (||) disregards the effect of the finite size of 
the vortex core, which removes the logarithmic infinity 
of B^{y) at r^, and thus reduces the amphtude of the 
higher Fourier components. At _B <C i3c2 this effect is 
accounted for by multipUcation of the London solution 
by a cutoff factor. Here a general remark seems 
appropriate. There is no general theory of B(r) valid 
at arbitrary temperature, and even if it existed (if the 
BCS-Gorkov-Eliashberg theory would apply and could 
be solved) the material parameters entering such a the- 
ory are not known with sufficient accuracy, e.g., the 
anisotropic electron mean free path /, the shape of the 
Fermi surface, and the coupling constant. Even when we 
use the Ginzburg-Landau (GL) theory to obtain a cutoff, 
we find that a general analytical solution does not exist, 
not even in the limit k — > c», which would be sufficient 
here. If the GL theory is applicable it applies down to 
B = {]. Below we derive the low-field cutoff factor from 
approximate analytical solutions of the GL theory and 
from a numerical solution. We are considering first an 
isotropic superconductor. 

The best analytical GL expression available was ob- 
tained by Clem for isotropic superconductors at low 
inductions B <C i?c2- Using a Lorentzian trial function 
for the order parameter |V'(r)P of an isolated vortex, 
Clem finds for large k ^ 1 



B,{G) 



^0 gKijg) 
S l + X^G"^' 
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V2aG^ + \-^Y'\ (3) 



Here Ki{x) = —Kq{x) is a modified Bessel function with 
the limits Ki{x) = 1/a;- (a;/2) ln(1.7139/x) (x < 1) and 
Ki{x) = (7r/2a;)i/2exp(-2;) (x > 1). From Eq. (|) we 
recover the London solution if the cores size shrinks to 
zero. The cutoff factor gKi{g) in Eq. (^) may be ap- 
proximated for all g values by exp(— \/2^G), or, less ac- 
curate but convenient for computations, by exp(— 2 ^^G^) 
as suggested in Ref. Q. The cutoff exp(— ^^G^/4) given 
in Ref. |jl^ was derived from the GL solution near Bc2, 
and is not valid at low B [B <^ ). At intermedi- 
ate fields the cutoff should interpolate between these two 
expressions. Therefore, the argument of the Gaussian 
cutoff used recently is smaller than the one we propose: 
1/4 111 or 1/2 instead of 2 vahd at low B. The cor- 
rect low-field cutoff yields a stronger field dependence of 
the SANS intensity than predicted for example in Ref. 

Clem's approximate analytical theory of the dilute vor- 
tex lattice was extended to larger fields and to anisotropic 
superconductors by Hao et al. |^ using the same type of 
variational approach. The resulting Fourier components 
for an isotropic superconductor may be written as 



B.{G) 



$0 



fooKi 



ifl + X^G^/' 



(4) 



where and foo are two variational parameters repre- 
senting the effective core radius of a vortex and the de- 
pression of the order parameter due to the overlap of vor- 
tex cores, respectively. For the cases of interest here (k > 
10) the two variational parameters have simple functional 
dependences on 6 = B/Bc2 and k : 



^„ = ^(V2 - ^) (1 + [1 _ 26(1 _ 6) 



211/2 



(5a) 



(5b) 



In Eqs. (|) <Po/S ^ B = bB^2 is the mean induction, 
which for 2bK^ > 1 may be equated to B^xt- 

For K ^ 1 the argument of Ki in the denominator of 
Eq. (|^) is much smaller than 1, thus we may use Ki(x) w 
1/x. Since for high Tc superconductors and typical i?oxt 
values, b is never larger than a few %, we may also neglect 
the field dependence of foo and , putting foo ~ 1 and S.v 
w v^^. For the analysis of measurements performed on 
heavy fermion superconductors, the field dependence of 
foo{b) can thus be disregarded (usually Bcxt < 1 T |^) 
but this may not be true for ^vib). For example, with 
UPts at Scxt = 1 T one has 6 « 0.4 (D and therefore 
« 0.854 X V2^. 

The smallest non-zero reciprocal vector for an equilat- 
eral triangular lattice is Gio = Gmin = o,l = (27r/S')a„ 
(see Fig. |l| for the definition of Uv), thus G^^jj^ = 
(87r2/V3)(B/$o)- This means that for the high Tc com- 
pounds at B ~ Boxt = 20 niT one has AG^^j^ ~ 10 > 1, if 
A^/^ = A = 1500 A is used. For UPta A is even larger [|). 
Accounting for the large value of AG^i„ = (47r/'\/3)6/v^ 
we may write 



B.iG) 



S AG2 



iCvG) Ki (e.G), 



(6) 



In this letter we test the applicability of formula (|6|) to 
recently published SANS results on CeRu2. 

The conventional superconductor CeRu2 has attracted 
some interest because of its complex phase diagram in the 
(_Boxt,r) plane. Notably, a reversible- irreversible line is 
observed. The form factor Bz{G) is easily obtained from 
the SANS cross-section [|o). The CeRu2 measurements 
of Bz{Gio) as a function of Bcxt are presented in Fig. ||. 
Because A is scalar, we derive from Eq. (^) 



^z(Gio) — 



31/4 VOT^^il^ 
A2 



27r\/2 

/27r%/2, , — 



t/$c 



(7) 



This expression depends only on the two parameters A 
and ^. The fits yield for the data recorded either in field 
cooling (FC) or zero field cooling (ZFC) procedure, A = 
1870 A and ^ = 84 A and A = 2090 A and ^ = 74 A, 
respectively. Taking the traditional point of view, the 
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FC data reflect the equilibrium properties of the vortex 
lattice. From these data k = 22 is larger than the pre- 
viously estimated n = 14.5 From the ^ value 
we compute Bc2 = <&o/(27r^^) = 4.7 T. Magnetization 
measurements at 1.8 K give Bc2 = 5.3 T |2^. The val- 
ues deduced from the FC neutron data are satisfactory 
in view of the well known difficulty to extract a reliable 
K value from magnetization measurements. 

The traditional Gaussian cutoff predicts ln[Bz{Gio)] 
cx -Bext, *-e- a straight line in Fig. |^. This is not observed. 

The generalization of Eq. to anisotropic penetra- 
tion length tensors reads for k >• 1 

Here uKi (u) is an anisotropic cutoff factor with 

u^ = 2{eGlHlGl){l + b^) \l-2bil-bf] , (8b) 



We shall not analyse the SANS data of UPta D with 
Eq. (||) because the conventional GL theory discussed 
here does not describe the phase diagram of this com- 
pound. We argue that the effect of the vortex cores in 
UPta is stronger than suggested by Joynt Q . 

We now consider the field distribution (probability) of 
the vortex lattice which is measured by fiSH p5[ and can 
be computed from the Fourier coefficients, Eq. (^). Its 
variance is = (i?^) — {Bz}'^, where (. . .) means the 
spatial average. One has 



2^ \B.{G)\^ 

G#0 



(12) 



Av separates into two factors, Av — Aq • fv{b) where 

$0 



Ao = 0.06092 



XxXy 



(13) 



is the London limit (^a;, 0) ||l7|] and /v(6) is a uni- 

versal function which accounts for the core size, 



uKi{u) w 1 - (uV4) ln(2.937/M2) for w < 1 . (8c) 



For the computation of Bz{G) we need to specify the 
geometry of the vortex lattice. As shown by Kogan [^ , 
for B ^ Bel, the angle characterizing this lattice (see 
Fig. [^) depends only on the penetration-length ratio : 

tana — V3{Xx/Xy). 



Using Kogan's formula (^), the form factor factorizes, 
Bz{Gpq) = Bo ■ bpq{b), where 



Bn^— — 



1 fsV^' $0 



TT^ \64/ XxXy 

and bpq{b) is a universal fmiction, 



(10) 



,{b)^{l-b^f-^f^^l^, (lla) 



1-26(1-6^ 



1/2 



1/2 



(lib) 



In Fig. 1^ we present 6io(6) computed from the varia- 
tional solution (11), the Gaussian cutoff (Ref. [19l) and 
the numerical solution of the GL equations ]23|] . Re- 
markably, the comparison between the variational and 
the numerical solutions shows that for b < 0.05 the first 
three Fourier coefficients Bz{G) deviate by < 10% and for 
6 < 0.01 by < 4%; even for b = 0.2 (0.3) the S^(Gio) (4) 
with (5) falls below the exact value by only 14% (18%), 
and even for small k — 5 this Clem-Hao approximation 
is reasonable. 



/2(6) = 0.12968 ' 

(p,g)#(0,0) 



pq ■ 



(14) 



cf. Fig. H The functions bio f^nd /v are very similar since 
in the sum ( p^ the six 6io equivalent terms dominate. 

Quite unexpectedly, the functions bio and /v are 
strongly field dependent even at low reduced fields 6, 
where the London model predicts constant bpq ~ 1 and 
(9) /v = 1- One has approximately \ — bpq(b) cx 1 — /v oc 6^/^, 



c/. Fig. ^. This finding is confirmed by the exact numer- 
ical solution of the GL theory |2^, depicted as dashed 
lines in Fig. |3|. This strong b dependence originates from 
the limit (|c) with = e^G^ = (87^/\/3)^>(G/Glo)^ 
which means that the cutoff factor uKi{u) is consider- 
ably less than unity except at very small b <C V^/ (Stt) — 
1/14.5 even for G = Gio. 

We are aware of only one investigation on a single crys- 
tal of the field dependence of the vortex lattice field dis- 
tribution ||]. From this ^SR study of YBa2Cu3 06.95 and 
the value 3^2 = 90 (10) T [||| we estimate Av « 5.04 
mT and 5.73 mT at B^^t = 0.5 T and 1.5 T, respectively. 
This leads to a ratio 7^oxp = Av(1.5 T)/Av(0.5 T) = 0.88 
while our computation (see Fig. ||) predicts T^-ql = 0.90. 
Therefore the conventional GL theory provides a simple 
and natural explanation of the observed field dependence 
of the observed field distribution in YBa2Cu3 06.95. 

In conclusion, we have shown that the effect of the 
finite core size on the Fourier components of the magnetic 
field in a conventional superconductor with large k is 
strong, even at low fields Bd < i? <C i?c2, since the 
cutoff factor in Eqs. (^ and (^) is uKi{u) < 1. This 
cutoff effect provides a natural explanation for recently 
published neutron and fiSK data without need to resort 
to unconventional theories. 
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FIG. 1. Definition of the primitive cell vectors &v and b„ 
and angle a of a distorted vortex lattice in real space, a is 
7r/2 minus the angle defined in Fig. 2 of Ref. Q. 
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FIG. 2. Form factor for the refiection [1,0] from the vortex 
lattice of CeRu2 as a function of the applied field. The points 
taken from Ref. have been obtained using either a field 
cooling or zero field cooling procedure. The lines are fits to 
Eq. (|). 



Figure 3, Yaouanc et al. 




0.0 



0.0 0.2 0.4 0.6 0.8 1.0 
b 



4 



FIG. 3. The universal functions 610(6) (Ha) (the largest 
reduced form factor, top) and /v(6) (14) (the reduced vari- 
ance, bottom) calculated in three ways: From this work (solid 
lines), from the Gaussian cutoff (dash-dotted linos), and from 
the exact Ginzburg-Landau solution (dashed lines). The in- 
serts plot these functions versus \fh to stretch the cusp-like 
h dependence of the correct cutoff at low reduced inductions 
6 = 0. Note the strong deviation of the previously used Gaus- 
sian from the correct cutoff. 



